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The isotopic evolution of the ground-state nuclear shapes and the systematics of one-quasiproton 
configurations are studied in neutron-rich odd-yl Yttrium and Niobium isotopes. We use a selfcon- 
sistent Hartree-Fock-Bogohubov formahsm based on the Gogny energy density functional with two 
parametrizations, DIS and DIM. The equal filling approximation is used to describe odd-A nuclei 
preserving both axial and time reversal symmetries. Shape-transition signatures are identified in 
the = 60 isotopes in both charge radii and spin-parities of the ground states. These signatures 
are a common characteristic for nuclei in the whole mass region. The nuclear deformation and shape 
coexistence inherent to this mass region are shown to play a relevant role in the understanding of the 
spectroscopic features of the ground and low-lying one-quasiproton states. Finally, a global picture 
of the neutron-rich A ~ 100 mass region from Krypton up to Molybdenum isotopes is illustrated 
with the systematics of the nuclear charge radii isotopic shifts. 

PACS numbers: 21.60.Jz, 21.10.Pc, 27.60.+j 
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I. INTRODUCTION 

The structural evolution as a function of the number 
of nucleons is a subject of increasing interest in nuclear 
structure, which is supported by a very intense activity 
on both theoretical and experimental sides (ll-fTsj. In 
particular, neutron-rich nuclei in the mass region around 
A — 100 have received special attention because of the 
interesting nuclear structure features merging there 0, 

Experimentally, the efforts are focused on different and 
complementary directions. The most relevant for the 
purpose of this work are related, on the one hand, to the 
mass determination 0-[lH in the case of exotic neutron- 
rich nuclei with the accuracy required for the modeling 
of astrophysical events [l^, [l^. In particular, the mass 
region of our concern here is highly significant to under- 
stand the nucleosynthesis path and the isotopic abun- 
dances generated by the astrophysical r process [l^l ■ On 
the other hand, the focus is on laser-spectroscopy exper- 
iments aimed to measure nuclear spins, magnetic dipole 
moments, spectroscopic quadrupole moments and mean- 
square charge radii from isotopic shifts. Considerable 
progress has been achieved in the last years (for a recent 
review, see Q and references therein), and special atten- 
tion has received the mass region studied in this work 

Theoretically, the region has been studied using 
phenomenological models [37l-l40| and microscopic ap- 
jroaches, based on the Relativistic Mean Field (RMF) 



41 [ well as nonrelativistic Skyrme [4^ |43| and Gogny 
Il-lil, H m energy density functional (EDF). All in 
all, the accumulated information on this mass region has 
allowed to establish some characteristic features, which 
can be associated to signatures of a shape transition at 
N = 60. 

Different nuclear properties sensitive to these struc- 
tural changes have been recently investigated [T^ - [l6| in 



several isotopic chains in this mass region. We used a 
selfconsistent Hartree-Fock-Bogoliubov (HFB) approxi- 
mation based on the finite range and density depen- 
dent Gogny-EDF [i^ and the equal filling approxima- 
tion (EFA) to deal with the odd nucleon. We analyzed 
bulk and spectroscopic properties of neutron-rich iso- 
topes with both even Z (Sr, Zr, and Mo isotopes) and 
odd Z (Rb isotopes). Our purpose in this paper is to 
complete the systematic study of the bulk and spectro- 
scopic properties in two chains of odd-Z isotopes, Yt- 
trium {Z = 39) and Niobium {Z = 41), as well as in the 
chain of Krypton isotopes {Z — 36). 

The description of odd- A nuclei involves additional dif- 
ficulties because the exact blocking procedure requires 
the breaking of time-reversal invariance, making the cal- 
culations more involved ^yrSlJ- In the present study we 
use the EFA, a prescription widely used in mean-field 
calculations to preserve the advantages of time-reversal 
invariance. The predictions arising from various treat- 
ments of the blocking have been studied in Ref. jso} , 
concluding that the EFA is sufficiently precise for most 
practical applications. More details of our procedure can 
be found in Ref. ^\^. 

In this work we consider two parametrizations of the 
Gogny-EDF, namely DIS [53 as the standard and most 



studied parameter set 



and DIM [EE] as 
the most recent effort to find a new parametrization that 
improves the predictions for masses maintaining the ex- 
cellent performance and predictive power of the former 
DIS. Our aim is first to verify the robustness of our pre- 
dictions with respect to the particular version of the EDF 
employed and second, to test the performance of DIM in 
the present context of the spectroscopy of odd-^ nuclei. 

The paper is organized as follows. In Sec. [TTl we 
present a brief description of the theoretical formalism 
used in the present work, i.e., the HFB-EFA framework. 
The results of our calculations for the considered nuclei 
are discussed in Sec. IIIIl where we pay attention to the 
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one-quasiparticle states and their spectroscopic evolution 
along the Y and Nb isotopic chains. We also compare our 
results with the available experimental data for charge 
radii and two-neutron separation energies. In Sec. IIVI 
we show the results for Krypton isotopes, as well as the 
systematics of the charge radii in the whole region from 
Krypton {Z = 36) up to Molybdenum {Z = 42). Finally, 
Sec. |V] is devoted to the concluding remarks and work 
perspectives. 

II. THEORETICAL FRAMEWORK 

Our theoretical framework to deal with odd-A nuclei 
is based on the Gogny-HFB-EFA formalism. In previous 
studies of even-even nuclei j6i, [Tj we have found advanta- 
geous to use the so called gradient method [s^l to obtain 
the solution of the HFB equations, leading to the (even 
number parity) vacuum |<I>). Within this method, the 
HFB equation is recast in terms of a minimization (vari- 
ational) process of the mean field energy. The Thouless 
parameters defining the most general HFB wave func- 
tions [ill are used as variational parameters. As it is cus- 
tomary in calculations with the Gogny force, the kinetic 
energy of the center of mass motion has been subtracted 
from the Routhian to be minimized in order to ensure 
that the center of mass is kept at rest. The exchange 
Coulomb energy was considered in the Slater approxi- 
mation and the contribution of the Coulomb interaction 
to the pairing field is neglected. Both axial and time- 
reversal are selfconsistent symmetries in our calculations 
for odd-A nuclei. Triaxial calculations have also been 
performed in the case of even-even Kr isotopes. 

The HFB ground-state wave function |<I>) of an even- 
even nucleus is defined by the condition of being the vac- 
uum of the annihilation quasiparticle operators of the 
Bogoliubov transformation [Sa, [s^. On the other hand, 
the ground and low-lying one-quasiparticle states of odd- 
A systems, like the ones considered in the present work, 
can be handled with blocked (odd number parity [ssLlsoj) 
HFB wave functions 

\^^B)=fiU\^). (1) 

where [Ib indicates the quasiparticle state to be blocked 
and stands for the indexes compatible with the symme- 
tries of the odd-nuclei, such as the angular momentum 
projection K and parity tt in the case of axial symmetry. 
As mentioned above, we use here the EFA to the ex- 
act blocking that preserves time-reversal invariance. In 
this approximation the unpaired nucleon is treated on 
an equal footing with its time-reversed state by sitting 
half a nucleon in a given orbital and the other half in 
its time-reversed partner. The microscopic justification 
has been first given in Ref. jiH using ideas of quantum 
statistical mechanics. The EFA energy can be obtained 
as the statistical average, with a given density matrix op- 
erator, and applying the variational principle to it, the 



HFB-EFA equation [i^ is obtained. The existence of a 
variational principle allows to use the gradient method to 
solve the HFB-EFA equation with the subsequent simpli- 
fication in the treatment of the constraints. 



The solution of the HFB-EFA (as well as the exact 
blocked HFB) equation depends upon the initial blocked 
level \xb ■ In the HFB-EFA case, the K quantum number 
is selfconsistently preserved along the calculation and so 
does the parity if octupole correlations are not allowed 
in the iterative process. Blocking levels with different 
K'" values lead to different quantum states of the odd-A 
nucleus, being the ground state the one with the lowest 
energy. One should notice that because of the selfconsis- 
tent nature of the whole procedure, for a given K'^ ^ there 
is no guarantee that the lowest energy solution for those 
K"^ values is obtained by taking the quasiparticle with 
the lowest energy as the initial blocked state. Therefore, 
several quasiparticles with the same K'^ must be con- 
sidered. In addition, in the present case and because of 
the presence of coexisting prolate, oblate, and spherical 
minima in some of the nuclei considered, blocked config- 
urations with those quadrupole deformations have to be 
explored. Constrained calculations have been performed 
to generate potential energy curves (PECs) for the even- 
even neighbor nuclei. One can find a systematic compi- 
lation of PECs obtained with Gogny-DlS in Ref. 0. In 
the case of odd-^ nuclei, the purpose of the computation 
of such PECs is twofold: First, they give us initial hints 
on the evolution of the different competing shapes in the 
considered nuclei and, second, they provide a whole set of 
prolate, spherical and oblate even-even HFB states (ref- 
erence states) for the subsequent treatment of the neigh- 
boring odd-A nuclei. In fact, once a reference (even- 
even) state with a given deformation is chosen, we use it 
to perform an additional (constrained) HFB calculation 
providing an unblocked fully-paired state corresponding 
to an odd average neutron number (false vacuum [47]) 
with the same deformation. Such prolate, spherical, and 
oblate false vacua are then used as input configurations 
in our subsequent blocking scheme (i.e., EFA). 



Thus, the minimization process has to be carried out 
several times, using different initial prolate, spherical and 
oblate (false) vacua. We have repeated each calculation, 
for a given false vacuum and K values from 1/2 up to 
15/2, using as initial blocking states the 12 quasiparti- 
cles corresponding to the lowest quasiparticle energies. 
The use of so many initial configurations is to guarantee 
that we are not missing the true ground state and all the 
lowest excited states. Note that for nuclei in this region 
of the nuclear chart, there exist several competing shapes 
at low excitation energy and therefore our procedure as- 
sures that the lowest energy solution can be reached for 
all values of the quadrupole moment Qiq and mass num- 
ber. 
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III. RESULTS FOR ODD-A YTTRIUM AND 
NIOBIUM ISOTOPES 



Being odd-Z nuclei, the spin and parity of odd-^ Yt- 
trium and Niobium isotopes are determined by the state 
occupied by the unpaired proton. The spectroscopic 
properties of the odd-^ isotopes are determined by the 
one-quasiproton configurations that, in principle, are ex- 
pected to be rather stable against variations in the num- 
ber of even neutrons. However, as it is known in neigh- 
boring nuclei, approaching TV ~ 60, the isotopes become 
well deformed [l|. 111] and the abrupt change in deforma- 
tion induces signatures in nuclear bulk properties like the 
two-neutron separation energies and the nuclear charge 
radii, as well as in spectroscopic properties. In particu- 
lar, the spin and parity of the nuclear ground state might 
flip suddenly from one isotope to another, reflecting the 
structural change. 



A. Low-lying one-quasiparticle states 

In Fig. [T]we can see the experimental excitation ener- 
gies and spin-parity assignments (a) in odd-^ neutron- 
rich Yttrium isotopes [61j . They are compared to 
the one-quasiproton states predicted by our Gogny-DlS 
HFB-EFA calculation (b), where the excited states in a 
given isotope are referred to the corresponding ground 
state, regardless its shape. Prolate configurations in our 
calculations are shown by black lines, oblate ones by red 
lines, and spherical ones by blue lines. The quasiparticle 
states are labeled by their K'^ quantum numbers. The 
most important configurations are joined by dashed lines 
following the isotopic evolution. In addition, the ground 
states are labeled by their asymptotic quantum numbers 
[N, n„A]K\ 

Experimentally, we observe J'^ = 1/2^ ground states 
in 87,89,91,93,95,97y^ then we have (5/2+) states in the 
heavier isotopes 99.ioi.io3,io5,io7y^ although their assign- 
ments are in these cases uncertain. The 9/2+ states ap- 
pear as excited states in the lighter isotopes, while 3/2~ 
and 5/2~ excited states are also found in most isotopes, 
but especially in the heavier ones. The most striking fea- 
ture observed is the abrupt transition at iV = 60 (A = 99) 
from 1/2^ to 5/2+ ground states. 

The theoretical interpretation of these findings can 
be understood from the analysis of our results in the 
lower figure (b). According to our calculations. Yttrium 
isotopes evolve from spherical shapes in 87-93y around 
N — 50, with the spherical pi/2 shell in the ground state 
and close 59/2 and /5/2 as excited states, to slightly de- 
formed oblate shapes in 95,97y^ g^^^^ finally to well de- 
formed prolate shapes in ^^-^OTy^ jj^ the lighter iso- 
topes the ground states correspond to 1/2^ states {pi/2), 
whereas the excited states correspond to the split lev- 
els coming from gg/2 and /5/2, with 9/2+ states as the 
lowest excited states in agreement with the measured 
low- lying spectra. The excited 9/2+ states decrease in 



energy as we move away from the magic neutron num- 
ber iV = 50 because of the incipient oblate deformation 
emerging, and they eventually become the ground state 
in ^''Y. In the case of the ^^Y isotope the oblate 9/2+ 
state is practically degenerate with the 1/2^ spherical 
state. In our description the nucleus ^^Y displays a 9/2+ 
oblate ground state with a l/2~ excited state at 0.7 MeV 
and an incipient prolate 5/2+ at 0.4 MeV that will be- 
come the ground state in the heavier isotopes ss-io^y 
All of these heavier prolate isotopes present oblate 9/2+ 
excited states at energies in the range 0.6 — 1.2 MeV. 
Thus, the experimentally observed leap of the ground- 
state spin-parity assignments at A'^ = 60 is well accounted 
for by the present calculations, which is interpreted as a 
sudden ground state shape change. 

Similarly to Fig. [U we show in Fig. [5] the correspond- 
ing results for Niobium isotopes. Experimentally 6l| we 
observe J'^ = 9/2+ ground states in the lighter ^s-sgjvjj^ 
isotopes with l/2~ excited states at close energy and 
then we find (5/2+) ground states in the heavier iso- 



topes 



101-105 



Nb, although as in the case of Y isotopes. 



the spin-parity assignments are still uncertain in these 
isotopes. A characteristic jump a.t N — 60 {A = 101) 
from 9/2+ to (5/2+) ground states is once more found. 
The calculations from Gogny-DlS (b) show in the lighter 
isotopes a clear correspondence between the measured 
ground state J'^ = 9/2+ and the calculated ones. The 
observed excited states l/2~ are also labeled in the cal- 
culations. They correspond to the excited states from 
the Pi/2 shell. However, at variance with experiment, 
our results indicate a transition in ^^Nb to the 7/2+ 
oblate state. This oblate configuration is kept all the 
way up to the heaviest isotope considered ^"^^Nb. The 
observed (5/2+) ground states are found in our calcu- 
lations as prolate configurations at excitation energies 
between 0.6 and 1 MeV in loi.ioa.iosjv^^^^ 

teresting to see whether the experimental ground states 
of ioi,i03,i05]\jj-, j^j^g confirmed to be 5/2+ states and to 
measure the heavier I0'''i09]\ji3^ where according to the 
calculations the 5/2+ states appear at a very high exci- 
tation energy relative to the ground state. 

The disagreement between the theoretical predictions 
and experimental data for the spin and parity of the heav- 
ier Nb isotopes could be understood if triaxiality effects 
are invoked. In Fig. [31 where the proton single particle 
levels for ^'^'^Zi are depicted as functions of the defor- 
mation parameter /3 [J], we observe for /3 = —0.2 (the 
position of the oblate minimum) a K"^ — 7/2+ orbital 
just above the Fermi level. The occupancy of this orbital 
produces the 7/2+ oblate ground state in ^"^Nb. The 
= 5/2+ orbital coming from the same 59/2 subshell 
and presumably responsible for the experimental spin 
and parity lies higher in energy. The situation is reversed 
in the prolate side where at /3 = 0.35 the = 5/2+ is 
below the = 7/2+ and is still above the Fermi level. 
The easiest way to connect the oblate and prolate sides is 
through the triaxial 7 degree of freedom, with 7 ranging 
from 60° (oblate side) to 0° in the prolate side. It is ob- 
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vious that the energy of the JC^ = 5/2+ orbital has to go 
below the energy of the = 7/2+ one at some 7 value 
and it is very likely that at that point the blocking of the 
= 5/2+ orbital will produce the ground state mini- 
mum. This is obviously a hand-waving argument as the 
K quantum number is not preserved along the 7 path but 
again, it is very likely, that both components K"^ = 5/2+ 
and K'^ = 7/2+ are going to be dominant in the orbital 
just above the Fermi level. The above argument calls for 
the necessity of carrying out a full triaxial calculation for 
the heaviest odd-Z Nb isotopes which is, for the moment, 
not possible as it will require an extension of our present 
computational codes for odd nuclei to include the role of 
triaxiality. Work along these lines is in progress. 

In the next figures we compare the spectroscopic prop- 
erties of the Gogny DIS and DIM parametrizations. In 
Figs. miSlini and[7]we have separated the prolate (a) and 
oblate (b) states and have plotted the most relevant hole 
states below zero energy and the particle states above. 
The ground states, either oblate or prolate, are indicated 
with a circle. Specifically, Figs. |4] and [5] correspond to Y 
isotopes with Gogny-DlS and Gogny-DIM, respectively, 
whereas Figs. |6] and [7] correspond to Nb isotopes with 
Gogny-DlS and Gogny-DIM, respectively. 

In Figs. 2] and [5] for Yttrium, in the prolate case (a) 
we have depicted the evolution of the 5/2+ ((79/2) states, 
which are ground states for ^^-^OTy^ ^j^e upper re- 
gion we find the 3/2-(/5/2), 5/2-(/5/2), and l/2+(d5/2), 
while in the lower region we have the 3/2+ ((/g/2), as it can 
be seen from Fig. [3] on the prolate side. Similar states 
are found in the prolate graphs (a) of Figs. |6] and [7] for 
Nb isotopes, but in this case the prolate configurations 
are not ground states. In the oblate case (b) we can see 
for Y isotopes the 9/2+(gg/2) state, which is ground state 
in ^^Y, the 7/2+(5g/2) as a particle state and the l/2~ 
and 3/2~ from as hole states, as it can be seen from 
Fig. [3] on the oblate side. The oblate configurations in 
Nb isotopes show the 7/2+(gg/2) as the ground state in 
all the isotopes depicted. In addition to the states shown 
for Y isotopes, we also show here the particle states 5/2+ 
and 3/2+ from gg/2- 

Very similar results are obtained from both 
parametrizations of the Gogny-EDF. The only dif- 
ference worth mentioning is that DIM produces slightly 
lower excited states and then a more compact level 
density. This feature can be understood from its larger 
effective mass that makes the single particle spectrum 
somewhat more dense with DIM. This answers our orig- 
inal question about robustness of the calculations and 
the reliability of DIM in what concerns the spectroscopic 
properties of the two parametrizations. 

To further illustrate the role of deformation and spin- 
parity assignments in the isotopic evolution, we display in 
Fig. |S]for Y (a) and Nb (b) isotopes, the axial quadrupole 
deformation /3 of the energy minima as a function of N 
with both DIS and DIM parametrizations. The defor- 
mation of the ground state for each isotope is encircled. 
We can see that beyond the semi-magic isotope with 



TV = 50 we start getting two minima in the prolate and 
oblate sectors. In the case of Y isotopes the spherical 
1/2^ states are ground states up to TV = 54 in agree- 
ment with experiment. For N — 52 and = 54 one 
can see the appearance of an oblate 9/2+ and a prolate 
1/2+ excited states, which are almost degenerate with 
the ground state. The next isotopes, with = 56 and 
N = 58, display a slightly oblate 9/2+ ground state and 
then suddenly a transition occurs at = 60 to well de- 
formed prolate 5/2+ ground states, in agreement with 
experiment. The calculations with both DIS and DIM 
are coincident. On the other hand, in the case of Nb iso- 
topes, we observe the spherical 9/2+ states being ground 
states along iV = 44 — 56 in agreement with experiment. 
However, starting at = 56 in the case of DIM and at 
A^ = 58 in the case of DIS, oblate 7/2+ ground states 
are obtained at variance with the experimental assign- 
ment (5/2+) for these states, which is nevertheless still 
uncertain, as already mentioned above. 

A qualitative understanding of the features just dis- 
cussed can be obtained by looking at the single particle 
levels for protons in ^°°Zr and depicted in Fig. [3] as a 
function of the beta deformation parameter. This Z = 40 
even-even nucleus is in between Y and Nb regarding the 
number of protons and therefore its single particle prop- 
erties should not differ too much (at least at a qualitative 
level) from the ones of Y and Nb. The insight obtained 
from Fig. [3] can only be qualitative owing to the selfcon- 
sistent character of the EFA that somehow incorporates 
the polarization effects induced by the single proton in 
the even-even core. 



B. Charge radii and two- neutron separation 
energies 

In Figs, ini Uni [TI] and [T^l we show the results for 
the charge radii differences (a) defined as (5(r^)^"'^ = 
(r^)^ — (r?)^°, calculated with the same corrections as 
in Ref . |14| , and for the two-neutron separation energies 
S2n (b) • Figures [3] and [TU] show the results for Yttrium 
isotopes calculated with Gogny-DlS and Gogny-DIM, 
respectively. Figures [TT] and [T^] are the corresponding 
ones for Niobium isotopes. Our results are compared 
with isotope shifts from laser spectroscopy experiments 
[13, m, '36] in the case of S{r'^) and with mass measure- 
ments from Refs. [Ull^l in the case of S2n- We plot the 
results corresponding to the spherical, oblate, and prolate 
shapes at the minima of the PECs. Results correspond- 
ing to the ground states are encircled in each isotope. 

The first thing to notice is the remarkable similarity be- 
tween the results obtained with Gogny-DlS and Gogny- 
DIM. The only visible difference between DIS and DIM 
is a somewhat better agreement with the experimental 
S'2,1 values in the case of DIM, as it can be expected 
from its improved fitting protocol, focusing in a more 
accurate reproduction of the nuclear masses 

For Y isotopes, the measured nuclear charge radii dif- 
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ferences depicted in (a) exhibit a sizable jump at A'' = 60 
where the radius suddenly increases. This observation is 
well accounted for in our calculations, where the encircled 
ground states show that a jump of the same magnitude 
occurs between = 58 and N — 60. As we can see, the 
origin of the jump is related to the shape transition from 
slightly oblate shapes to well deformed prolate configu- 
rations beyond A^ = 60. The charge radius demonstrates 
once more to be a nuclear property very sensitive to those 
shape variations. It is also worth mentioning that the 
sudden change of the radius is perfectly correlated with 
the change in the spin-parity of the ground states studied 
above. 



In the case of Nb isotopes, shown in Figs. [TT] and [T^l 
the data are still scarce [33], but nevertheless a jump 
between A^ = 58 and A^ = 60 is also experimentally ob- 
served, although not as sharp as in the previous case. 
This is also associated with the transition from 9/2+ to 
(5/2+) observed in these isotopes between A^ = 58 and 
A^ = 60. As we have mentioned, the calculations produce 
only a gradual transition from spherical to oblate shapes, 
but the prolate shapes associated to the 5/2+ states never 
become ground states. In order to follow the experimen- 
tal trend we will need a transition at A^ = 60 to get our 
calculations to the data points, which are explained in 
any case by the prolate shapes. We faced in the past (T^ | 
similar problems in the case of Mo {Z — 42) isotopes, 
which are the neighbors of Nb {Z = 41). In that case 
it was demonstrated that, for even-even Mo isotopes, 
the incipient emergence of triaxiality was at the origin 
of the experimentally observed radius evolution. The 
quadrupole deformations of the triaxial minima were in 
Mo isotopes very close to the prolate ones and therefore, 
their effect on the radii where also similar to the prolate 
case. Similar type of arguments could be use here for 
Nb isotopes and one would expect the triaxial degrees of 
freedom to be more involved, but our present technical 
capability prevents us to carry out triaxial calculations 
for odd-A nuclei. 



The results for S2n agree in general with the measure- 
ments, but we get systematically a shell effect at A^ = 50 
stronger than experiment. This is a well known feature 
of any mean field approach 1435 j which is cured once be- 
yond mean-field configuration mixing calculations, in the 
spirit of the Generator Coordinate Method (GCM), are 
implemented in the formalism. In any case, it would be 
very helpful to extend and improve mass measurements, 
reducing the still large uncertainties in neutron-rich iso- 
topes. Most of the reported measurements in this ex- 
otic mass region are based on /3-endpoint measurements, 
which are not completely reliable because they lead to 
very strong binding |3l| . 



IV. CHARGE RADII SYSTEMATICS IN THE 
KR-MO REGION 



In this section we discuss the isotopic evolution of the 
charge radii in the whole region from Kr (Z = 36) up to 
Mo {Z — 42). These isotope shifts have been found to be 
very sensitive probes of nuclear shape transitions and it 
is worth studying globally the whole Kr-Mo region dis- 
cussing the similarities and differences among the various 
isotopic chains. 

In Fig. [13] we have compiled the measured isotope 
shifts 5(rj) for Mo Nb 0, Zr [1^, Y 0, HI], Sr 
[l^*], Rb [21|, and Kr 23] isotopes. These data are shown 
with solid symbols and connected with continuous lines. 
In addition, we can also find our results from Gogny- 
DIS-HFB calculations, which are shown with open sym- 
bols connected with dashed lines. Very similar results 
are provided by the Gogny-DIM EDF. As we can see in 
Fig. 1131 this mass region is characterized by a jump of 
the mean square charge radii at around A^ = 60, which 
refiects the sudden increase of deformation or shape tran- 
sition that occurs for these isotones. The increase in the 
charge radii is maximum for Y isotopes in the middle re- 
gion and it corresponds to a change from oblate to pro- 
late shapes at A^ = 60, as it was shown in the previous 
section. Around Y nuclei, we find that Zr and Sr iso- 
topes also show this big effect that was studied in Refs. 
(14 [T^ and interpreted as an oblate-prolate transition as 
well. The same is also true for Rb isotopes, as discussed 
in Ref. Ilal. T he behavior of Mo isotopes was studied in 
Refs. [Tj, [13, where it was shown that the suppression 
of the jump has its origin in the onset of triaxiality that 
smooths the evolution of the radii. The case of Nb iso- 
topes studied in this paper shows strong similarities with 
the case of Mo isotopes and triaxiality is likely again the 
reason for the discrepancies. 

Finally, we also include results for Kr isotopes, mea- 
sured in Ref. [l^- The isotope shifts for Kr isotopes do 
not show any abrupt change, but increase smoothly with 
the neutron number N. To understand this behavior in Kr 
nuclei, at variance with the abrupt change observed in the 
heavier neighboring isotopic chains, we have studied the 
evolution of the corresponding Potential Energy Surfaces 
(PESs) (i.e., Q — J energy contour plots) in even-even Kr 
isotopes with the help of constrained HFB calculations 
along the lines discussed in Refs. [1, [t*!. We find that the 
lighter Kr isotopes present very shallow PESs centered at 
the spherical shapes. Heavier Kr nuclei become gradually 
oblate in their ground states with prolate configurations 
at higher energies. Therefore, the oblate configurations 
are stabilized in Kr isotopes and we do not find any tran- 
sition to prolate shapes. The result is a smooth variation 
of the ground-state structure and, as a consequence, of 
the isotope shifts. To illustrate this point we show in 
Fig. [14] the triaxial landscapes for even-even Kr isotopes 
in the vicinity of A^ = 60. We can see that the ground 



state in ^ Kr (a) is oblate with a very soft variation in Q 
and especially in the 7 direction. The nucleus ^"^Kr (b) 
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develops an oblate ground state, which is again 7-soft. 
A prolate saddle point is also apparent. The heavier iso- 
topes ^^Kr (N = 60) (c) and ^^Kr {N = 62) (d) develop 
two well defined oblate and prolate minima, separated by 
energy barriers in the Q and 7 degrees of freedom. The 
ground state corresponds always to the oblate shapes, 
while the prolate shapes appear about 1 MeV higher in 
energy. 

In general, our theoretical calculations describe suc- 
cessfully this phenomenology. The smooth behavior ob- 
served in Kr isotopes is a consequence of the stabilization 
of the oblate shapes along the isotopic chain. The sud- 
den increase of the charge radii at TV = 60 observed in 
Rb, Sr, Y, and Zr isotopes is explained by the oblate- 
to-prolate shape transition. Finally, the tendency to a 
smooth behavior observed again in Nb and Mo isotopes 
is explained by the onset of a region of triaxiality. 



V. CONCLUSIONS 

In this work we have studied the shape evolution in 
odd-A Yttrium and Niobium isotopes from microscopic 
selfconsistent Gogny-EDF HFB-EFA calculations. We 
have also analyzed the isotopic evolution of the one- 
quasiproton configurations, comparing the predictions of 
the parametrizations DIS and DIM of the Gogny-EDF 
and demonstrating the robustness of our calculations. 

Two-neutron separation energies, charge radii, and the 
spin-parity of the ground states have been analyzed in a 
search for signatures of shape transitions. We have found 
that the charge radii and the spin-parity of the ground 
states are very sensitive to shape changes in these iso- 
topes. In addition, the signatures found are all consis- 
tent to each other. We have also shown that the qual- 
ity of the spectroscopic results obtained with the recent 
incarnation DIM of the Gogny-EDF is comparable to 
the one obtained with the standard parametrization DIS. 
We conclude that both DIM and DIS parametrizations 
reproduce, at least qualitatively, the main features ob- 
served in the isotopic trends of the considered neutron- 
rich isotopes. 

The observed spin-parity of the ground states exhibits 
a sudden change between TV = 58 and = 60 isotopes in 
both isotopic chains. This is a clear signature of a nuclear 
structure transition, which is correlated to the observa- 
tion of a similar change in the charge radii. According 
to our calculations this structural change is explained by 
a shape transition taking place at this neutron number. 
The relevant feature to stress is that the jump in both 
spin-parity and charge radius is a signature of a shape 
transition, where the systematics is suddenly broken to 
a new situation. In the case of Y isotopes, the exper- 
imental spin-parities are well reproduced, except in the 
case of A^ = 58, where the calculation predicts a transi- 
tion going across an oblate configuration, apparently not 
seen experimentally. In the case of Nb isotopes, the pre- 
dicted transition is to oblate states, while the experimen- 



tal spin-parities of heavier isotopes do not support these 
results. Since this is a region of an emergent triaxial- 
ity, as it has been shown in neighboring even-even nuclei, 
this is a plausible explanation of the discrepancy. Fi- 
nally, one should also take present that the experimental 
spin-parity assignments in the heavier isotopes of both 
chains are either uncertain or taken from systematics, 
and therefore, subject to revision. 

We have performed a systematic study of the behav- 
ior of the charge radii in the region of neutron-rich iso- 
topes from Kr {Z — 36) up to Mo {Z — 42) to stress 
the manifest enhanced sensitivity of these observables to 
shape transitions. Comparison with the available data 
from laser spectroscopy demonstrates the quality and ro- 
bustness of the Gogny-HFB description that is able to 
reproduce the main features of this behavior and offers 
an intuitive interpretation in terms of sharp or soft shape 
transitions, as well as triaxiality. 

The experimental information available in neutron-rich 
Y and Nb isotopes, and in general in this mass region, is 
still incomplete. It would be desirable to extend the ex- 
perimental programs for masses, charge radii, and spec- 
troscopic measurements to these exotic regions at exist- 
ing facilities like the tandem Penning trap mass spec- 
trometer ISOLTRAP [III at the On-Line Isotope Mass 
Separator (ISOLDE) facility at CERN and the Ion Guide 
Isotope Separator On-Line (IGISOL) facility [ll| at the 
University of Jyvaskyla, or at future ones like the preci- 
sion measurements of very short-lived nuclei using an ad- 
vanced trapping system for highly- charged ions(MATS) 
and laser spectroscopy (LaSpec) at the Facility for 
Antiproton and Ion Research (FAIR) , where we can learn 
much about structural evolution in nuclear systems. 

Although the present theoretical approach explains 
reasonably well the basic features observed in the data, 
theoretical efforts should be still pushed forward by im- 
proving the formalism by including triaxial degrees of 
freedom in odd- A isotopes, in particular in those regions 
where this can be an issue, or by dealing with the odd sys- 
tems with exact blocking treatments. The quality of our 
mean field description could be also improved from con- 
figuration mixing calculations in the spirit of the GCM 
with the quadrupolc moment as collective coordinate. It 
is already known 43] that such a configuration mixing 
reduces the jump in the predicted S2n, as compared to 
pure mean-field predictions, when crossing shell closures, 
thus improving the agreement with experiment. This 
could be particularly relevant for the light isotopes con- 
sidered in the present study, where the spherical minima 
are rather shallow. For heavier isotopes, the two minima, 
oblate and prolate, are separated by spherical barriers of 
about 3 MeV and appear typically about 1 MeV apart. 
The effect here is not expected to be significant because 
a single shape would be enough to account for the prop- 
erties studied. 
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FIG. 1: (Color online) Experimental (a) excitation energies and spin-parity assignments of the non-collective states are 
compared with HFB-EFA results (b) for the one-quasiproton states in odd-j4 Yttrium isotopes. Prolate configurations are 
shown by black lines, oblate ones by red lines, and spherical ones by blue lines. 



89 91 93 95 97 99 101 103 105 107 109 

(a) 
Nb - exp 



(7/2 ) 7/2 ■ 



3/2 7/2 



+ 7/2 



(1/2") HI- 1/2" 



3/2 

(5/2^) r-(7/2^5/2") 

t;7 (5^ 

(3/2) (l^i") 



(9/2^) 9/2^ 9/2^ 9/2^ 9/2^ 9/2^ (5/2^) (5/2^) (5/2^) 



1/2" 



1/2" 



1/2^ 



3/2^'' 3/2" 



1/2^ - 



3/2_ 3/2" 



3/2" / ' 
'1/2 



,'5/2" 



-5/2^ 



3/2" .^^Z;^ 3/2" 



3/2 



+ / 3/21,' 3/2" 



5/2 /,? 3/2" 



1/2^ 



9/2 



J 'I 

-a 9/2" 



9/2 



^ ' 1/2" 



i5/2 
3/2" 



1/2 ; 

5/2' 



5/2"'" 



.'if 1/2 



1/2 



5/21: 
3' 



; 5/2" 



.,^5/2" '5/2" 



^ 3/2 



1/2" '\ 3/2"; 



\ \ /II 

I II 



5/2^ 



(b) 



Nb - DIS 



7/2_ 7/2_ 7/2_ 7/2_ 7/2_ 7/2 

^ 9/2^ 9/2^ 9/2^ 9/2"[404]7/2^ [404]7/2^ [404]7/2^ 

9/2 §9/2 §9/2 §9/2 §9/2 



9/2 



89 91 93 



'=9/2 
95 



97 



[404J7/2 [404J7/2 [404J7/2 
99 101 103 105 107 109 



A 



FIG. 2: (Color online) Same as in Fig. [Jl but for Niobium isotopes. 
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FIG. 3: (Color online) Proton single particle energies for ^""Zr as a function of the quadrupole deformation parameter /?. 
The Fermi level is plotted by a thick dotted line. The results correspond to the Gogny-DlS EDF. Solid lines correspond to 
levels with positive parity, whereas dashed lines correspond to negative parity states. Asymptotic Nilsson quantum numbers 
[N,nz, A.]K'^ are also shown at the vertical lines where the minima of the potential energy curves are located. 
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FIG. 5: (Color online) Same as in Fig. H but for Gogny-DIM. 
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FIG. 6: (Color online) Same as in Fig. IH but for Niobium isotopes. 
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FIG. 7: (Color online) Same as in Fig. [5l but for Niobium isotopes. 
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FIG. 8: (Color online) Isotopic evolution of the quadrupole deformation parameter /3 of the energy minima obtained from 
Gogny DIS and DIM calculations for Yttrium (a) and Niobium (b) isotopes. 
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FIG. 10: (Color online) Same as in Fig. |9l but for Gogny-DIM. 




11: (Color online) Same as in Fig. |9l but for Niobium isotopes. Experimental data for radii are from Ref. 




FIG. 12: (Color online) Same as in Fig. |9] but for Niobium isotopes and Gogny-DIM. 
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FIG. 13: (Color online) Gogny-DlS HFB results for 5{r^) compared to the measured values for Kr, Rb, Sr, Y, Zr, Nb, and Mo 
isotopic chains. 
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FIG. 14: (Color online) Q - 7 planes for ^^Kt (a), ^''Kr (b), ^®Kr (c), and ^*Kr (d) with the Gogny-DlS EDF. The absolute 
minimum is marked by a bullet. The full contour lines correspond to energies ec (relatives to the absolute minimum) of 0.25, 
0.75, 1.25, 1.75, 2.25 and 2.75 MeV. Close to those full line contours other contour lines corresponding to ec + 0.1 MeV are 
also depicted. The purpose of these dotted contour lines is to give the direction of increasing energy as well as a visual idea of 
the corresponding slope. Finally, the dashed contour lines corresponding to energies from 4 MeV up to 10 MeV in steps of 1 
MeV are depicted to mark the region where the potential energy starts to grow rapidly. 



